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Abstract 



> , 

, We analyze the dynamics of dissipation and relaxation in the unbroken and 



broken symmetry phases of scalar theory in the nonlinear regime for large 
initial energy densities, and after linear unstabilities (parametric or spinodal) 
are shut-off by the quantum backreaction. A new time scale emerges that sep- 
. arates the linear from the non-linear regimes. This scale is non-perturbative in 



the coupling and initial amplitude. The non-perturbative evolution is stud- 
0^1 led within the context of the 0{N) vector model in the large limit. A 

Q^' combination of numerical analysis and the implementation of a dynamical 

^ ■ renormalization group resummation via multitime scale analysis reveals the 

presence of unstable bands in the nonlinear regime. These are associated 
with power law growth of quantum fluctuations, that result in power law 
^ \ relaxation and dissipation with non-universal and non-perturbative dy- 

' namical anomalous exponents. We find that there is substantial particle 

production during this non-linear evolution which is of the same order as 
that in the linear regime and results in a non-perturbative distribution. The 
expectation value of the scalar field vanishes asymptotically transferring all 
of the initial energy into produced particles via the non-linear resonances 
in the unbroken symmetry phase. The effective mass squared for the quan- 
tum modes tends asymptotically to a constant plus oscillating 0{\/t) terms. 
This slow approach to asymptotia causes the power behaviour of the modes 
which become free harmonic modes for late enough time. We derive a simple 
expression for the equation of state for the fluid of produced particles that 
interpolates between radiation-type and dust-type equations according to the 
initial value of the order parameter for unbroken symmetry. For broken sym- 
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metry the produced particles are asymptotically massless Goldstone bosons 
with an ultrarelativistic equation of state. We find the onset of a novel form 
of dynamical Bose condensation in the collisionless regime in the absence of 
thermalization. 



I. INTRODUCTION AND MOTIVATION 

The next generation of high luminosity heavy ion colliders at Brookhaven and CERN will 
offer the possibility of probing the dynamics of states of high energy density and possibly 
strongly out of equilibrium. The energy densities attained for central collisions at central 
rapidity will hopefully allow to study the quark-gluon plasma and also the chiral phase 
transition in a situation that parallels that achieved in the very early stages of the Universe 
0" [H- Dynamical phenomena, non-equilibrium and collective effects are expected to take 
place on time scales of a few tens of fm/c and length scales of few fermis. 

This unparalleled short time and length scale regime for dynamical phenomena soon 
to be probed experimentally, has sparked a considerable effort to study the dynamics of 
strongly out of equilibrium situations within the realm of quantum field theory. 

The usual semi-phenomenological framework to study the dynamics is based on the 
transport approach in terms of single (quasi) particle distribution functions with collisional 
relaxation [§. The best known dynamical processes of relaxation are those of (few body) 
collisions and dephasing processes akin to Landau damping [^- |T^-|T6[. Our understand- 



ing of these relaxational processes is usually based on perturbative expansions, linearized 
approximations or small departures from equilibrium. The validity of these coarse grained 
descriptions of relaxational dynamics within the realm of high energy and high density 
regimes in quantum field theory is not clear and a closer scrutiny of relaxational phenomena 
is warranted. 

Whereas equilibrium phenomena are fairly well understood and there are a variety of tools 
to study perturbative and non-perturbative aspects, strongly out of equilibrium phenomena 
are not well understood and require different techniques. 

Our goal is to deal with the out of equilibrium evolution for large energy densities in 
field theory. That is, a large number of particles per volume m~^, where m is the typical 
mass scale in the theory. The most familiar techniques of field theory, based on the S-matrix 
formulation of transition amplitudes and perturbation theory apply in the opposite limit of 
low energy density and since they only provide information on in out matrix elements, 
are unsuitable for calculations of time dependent expectation values. 

Recently non-perturbative approaches to study particle production |T^, dynamics of 
phase transitions [p!7| , p!8| and novel forms of dissipation ||20| have emerged that provide a 



promising framework to study the dynamics for large energy densities like in heavy ion 
collisions. 

Similar tools are also necessary to describe consistently the dynamical processes in the 
Early Universe [Q. In particular it has been recognized that novel phenomena associated 
with parametric amplification of quantum fluctuations can play an important role in the 
process of reheating and thermalization [|T^,^,^. It must be noticed that the dynamics in 
cosmological spacetimes is dramatically different to the dynamics in Minkowski spacetime. 
Both in fixed FRW and de Sitter backgrounds and in a dynamical geometry ||^ 
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the dynamical evolution is qualitatively and quantitatively different to the Minkowski case 
considered in the present paper. 

Our program to study the dynamical aspects of relaxation out of equilibrium both in 
the linear and non-linear regime has revealed new features of relaxation in the collisionless 
regime in scalar field theories [^,0 • Recent investigation of scalar field theories in the non- 
linear regime, including self-consistently the effects of quantum backreaction in an energy 
conserving and renormalizable framework have pointed out to a wealth of interesting non- 



perturbative phenomena both in the broken and unbroken symmetry phases |2^. These 
new phenomena are a consequence of the non-equilibrium evolution of an initial state of 
large energy density which results in copious particle production leading to non-thermal 
and non-perturbative distribution of particles. Our studies have focused on the situation 
in which the amplitude of the expectation value of the scalar field is non-perturbatively 
large, A ^ a/A < $ > /m ^ ("^ is the mass of the scalar field and A the self- 

coupling) and most of the energy of the initial state is stored in the 'zero mode', i.e. the 
(translational invariant) expectation value of the scalar field $. Under these circumstances 
the initial energy density e ~ rn^/X. During the dynamical evolution the energy initially 
stored in one (or few) modes of the field is transferred to other modes resulting in copious 
particle production initially either by parametric amplification of quantum fluctuations in 
the unbroken symmetry phase, or spinodal instabilities in the broken symmetry phase. This 
mechanism of energy dissipation and particle production results in a number of produced 
particles per unit volume M oc m?/\ which for weak coupling is non-perturbatively large 
||2CI|| . We call 'linear regime' to this first stage dominated either by parametric or spinodal 
unstabilities. 



We recognized a new dynamical time scale ti where the linear regime ends. By the 
time ti the effects of the quantum fluctuation on the dynamical evolution become of the 
same order as the classical contribution given by the evolution of the expectation value of 
the field. The 'non-linear' regime starts by the time ti. In the case of broken symmetry, 
this time scale corresponds to the spinodal scale at which the backreaction of quantum 
fluctuations shut-off the spinodal instabilities. At this scale non-perturbative physics sets 
in and the non-linearities of the full quantum theory determine the evolution. This time 
scale ti, which we call the non-linear time, is a non-universal feature of the dynamics and 
depends strongly on the initial state and non-perturbatively on the coupling, as ti oc log[A~^] 
for weak coupling pO[| . 

The purpose of this paper is to carefully analyze the nonlinear dynamics of relaxation 
after the time ti in a weakly coupled scalar field theory within a non-perturbative self- 
consistent scheme. We focus on the asymptotic time regime both in the unbroken and 
broken symmetry states. Here we provide refined numerical analysis of the non-equilibrium 
evolution that reveal the onset of widely separated relaxational time scales. We use a 
dynamical renormalization group implemented via a multitime scale analysis to provide an 
analytic description of the asymptotic dynamics and establish that relaxation occurs via 
power laws with anomalous dynamical exponents. 

The main results of this article can be summarized as follows: 

• The hierarchy of separated time scales allow us to implement a dynamical renormal- 
ization group resummation via the method of multitime scale analysis. The novel 
result that emerges from this combination of numerical and dynamical renormaliza- 
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tion group analysis is the presence of non-linear resonances that lead to asymptotic 

relaxation described by non-universal power laws. These power laws are determined 
by dynamical anomalous exponents which depend non-perturbatively on the coupling. 

• The effective mass felt by the quantum field modes in high energy density situations 
varies with time and depends on the fields themselves reflecting the nonlinear character 
of the dynamics. Both for broken and unbroken symmetry the effective mass tends 
asymptotically to a constant. This constant is non-zero and depends in the initial 
state for unbroken symmetry. For broken symmetry, the effective mass tends to zero 
corresponding to Goldstone bosons. In both cases the effective mass approaches its 
t — oo value as 1/t times oscillating functions. 

• The fact that the effective mass tends asymptotically to a constant imphes that the 
modes becomes effectively free. Non-resonant modes oscillate harmonically for times 
t > ti. Resonant modes change from non- universal power behaviour to oscillatory 
behaviour at a time that depends on the wavenumber of the mode. Only the modes 
in the borders of the band resonate indefinitely. 

• In the unbroken symmetry case, we find that the expectation value relaxes to zero 
asymptotically with a non-universal power law. The initial energy density which is 
non-perturbatively large goes completely into production of massive particles. The 
asymptotic particle distribution is localized within a band determined by the initial 
conditions with non-perturbatively large amplitude ~ 1 /y/X which could be described 
as a 'semiclassical condensate' in the unbroken phase. 

• The particle distribution in the condensate is nontrivial. We establish sum rules that 
yield explicit values for integrals over such asymptotic distribution. We derive in this 
way the asymptotic equation of state. For unbroken symmetry, it interpolates between 
dust and radiation according to the initial field amplitude. 

• In the broken symmetry phase when the initial expectation value of the scalar field 
(order parameter) is close to the false vacuum (at the origin) we similarly find the 
onset of non-linear resonances at times larger than the non-linear (spinodal) time ti. 
The expectation value of the order parameter approaches for late times a nonzero limit 
that depends on the initial conditions. 

• The effective time dependent mass vanishes as 1/t resulting in that the asymptotic 
states are Goldstone bosons. We also find a hierarchy of time scales of which ti oc 
ln[A^^] is the first an another longer time scale t2 oc 1/\/A. As a consequence of the 
non-linear resonances the particle distribution becomes localized for t > ti at very 
low momentum resulting again in a 'semiclassical condensate' with non-perturbatively 
large amplitude ~ 1/\/A. Asymptotically the equation of state is that of radiation 
although the particle distribution (Goldstone bosons) is non-thermal. 

• For even larger time scales, t ~ \/V (where V stands for the volume of the system), we 
find for broken symmetry a novel form of Bose condensation in the coUisionless regime 
that results from a linear growth in time of homogeneous quantum fiuctuations. 
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The article is organized as follows: in section II we briefly summarize the nature of the 
approximations, the non-equilibrium framework and some of the previous results for the 
beneflt of the reader and for coherence. In section III we study the unbroken symmetry 
case and distinguish the linear regime of parametric ampliflcation t <ti from the non-linear 
regime {t > ti) in which the backreaction of quantum fluctuations dominates the evolution. 
In section IV we study the latter regime in the unbroken symmetry case and we flnd that 
particle production continues beyond t > ti and non-linear resonances develop leading to 
power law relaxation. We provide a full numerical analysis and implement a renormalization 
group resummation of secular terms via a multi-timescale approach. Asymptotic sum rules 
and the equation of state are discussed in detail. In section V we study the dynamics in the 
broken symmetry phase establishing a difference between the early and intermediate scales 
dominated by spinodal instabilities and the asymptotically large time scale dominated by 
non-linear resonances leading to power law relaxation. We provide a numerical analysis 
as well as arguments based on multitime scale resummation. We flnd a novel form of 
Bose condensation with a quadratic time dependence for the formation of an homogeneous 
condensate. 

Conclusions and further questions are summarized at the end of the article. 



II. PRELIMINARIES 



As our previous studies of scalar fleld theory have revealed ||2^, there are two very 
important parameters that influence the quantum dynamics: the strength of the coupling 
constant A and the initial energy density in units of the scalar fleld mass m. If in the initial 
state most of the energy is stored in few modes, the energy density is determined by the 
amplitude of the expectation value of those modes A \/X < $ > /m. The value of 
this fleld amplitude determines the regime of applicability of perturbation theory methods. 
Usual S-matrix theory treatment (in terms of a perturbative expansion) are valid in the 
small amplitude regime A «1 even for high energies. However, when the initial state has 
a large energy density perturbative methods are invalid. 

We shall be concerned here with the non-perturbative regime in which A = a/A < $ > 
/m ~ 0{1). It is important to point out that for large fleld amplitude, even for very weakly 
coupled theories non-linear effects will be important and must be treated non-perturbatively. 
This is the case under consideration. Having recognized the non-perturbative nature of the 
problem for large amplitudes we must invoke a non-perturbative, consistent calculational 
scheme which respects the symmetries (continuous global symmetries and energy-momentum 
conservation), is renormalizable and lends itself to a numerical treatment. 



We are thus led to consider the 0{N) vector model with quartic self- interaction pO[ and 
the scalar fleld in the vector representation of 0{N). 
The action and Lagrangian density are given by, 

S = I d'^x £, 
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d^^{x) - V{^{x)) 
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The canonical momentum conjugate to is, 

n(x) = ${x), (2.2) 

and the Hamiltonian is given by, 

H{t) = J d'x[ln\x) + \ [V$(x)]2 + . (2.3) 



The calculation of expectation values requires the study of a density matrix whether 
or not the initial state is pure or mixed. Its time evolution in the Schrodinger picture is 
determined by the quantum Liouville equation 

^f = [^,P]- (2.4) 
The expectation value of any physical magnitude A is given as usual by 

< A > (t) = Ti[p{t) A] . (2.5) 



The time evolution of all physical magnitudes is unitary as we see from eq. 

In the present case we will restrict ourselves to a translationally invariant situation, i.e. 
the density matrix commutes with the total momentum operator. In this case the order 
parameter < $(x, t) > will be independent of the spatial coordinates x and only depends 
on time. 

We write the field $ as $ = (cr, vf) where vr represents the — 1 'pions', and choose 
the coupling A to remain fixed in the large limit. In what follows, we will consider two 
different cases of the potential ( |2.1| ) V{a,TT), with (m^ < 0) or without (m^ > 0) symmetry 
breaking. 

We can decompose the field a into its expectation value and fluctuations xi^^t) about 

it: 



a{x,t) = (f){t) VN + x{x,t) , (2.6) 

with (f){t) being a c-number of order one in the N —>■ oo limit and x an operator. 

To leading order the large A^-limit is implemented by considering a Hartree-like factor- 
ization (neglecting 1/A^terms) and assuming 0{N — 1) invariance by writing 

N-l 



n{x,t) = tp{x,t) {1,1, ■■■,!) (2.7) 
where il){x,t) is a quantum operator pO|. Alternatively the large A^ expansion is systemati- 



cally implemented by introducing an auxiliary field |17|. To leading order the two methods 
are equivalent. 

The generating functional of real time non-equilibrium Green's functions can be written 
in terms of a path integral along a complex contour in time, corresponding to forward 
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and backward time evolution and if the initial density matrix describes a state of local 
thermodynamic equilibrium at finite temperature a branch down the imaginary time axis. 
This requires doubling the number of fields which now carry a label ± corresponding to 
forward (+), and backward (— ) time evolution ||25| . |26[ . 

We shall not rederive here the field evolution equations for translationally invariant 
quantum states the reader is referred to the literature. (See refs. p5|-p7|). In the leading 
order in the large N approximation the theory becomes Gaussian at the expense of a self- 
consistent condition |[T7| . pO| , |27[| , this in turn entails that the Heisenberg field operator ip{x,t) 
can be written as 



d^k 1 
(2^71 



(2.8) 



where , a\ are the canonical creation and annihilation operators, the mode functions fk{t) 



are solutions of the Heisenberg equations of motion [ll7| , |20| , ^ to be specified below for each 
case. 

Our choice of initial conditions on the density matrix is that of the vacuum for the 
instantaneous modes of the Hamiltonian at the initial time [^,0 • Therefore we choose the 
initial conditions on the mode functions to represent positive energy particle states of the 
instantaneous Hamiltonian at t 



0, which is taken to be the initial time. That is, 



A(o) 



A(0) = -l^ 



(2.9) 



where the mass Mq determines the frequencies uJk{0) and will be defined explicitly later as 
a function of 0(0) [see eqs.( p.5|) and (|5.4| )]. 

With these boundary conditions, the mode functions /fc(0) correspond to positive fre- 
quency modes (particles) of the instantaneous quadratic Hamiltonian for oscillators of mass 
Mo. 

We point out that the behaviour of the system depends mildly on the initial conditions 
on the mode functions as we have found by varying eqs. (|2.9| ) within a wide range. In 
particular, the various types of linear and nonlinear resonances are independent of these 



initial conditions 20,27 



It proves convenient to introduce the following dimensionless quantities: 

__k_ r. _Wj^ 

\m\ \m\ 

A 



r 
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A 



2|m|2 
A 

2|mP 



m\ fkit) . 



S(0) = 



(2.10) 

(2.11) 

(2.12) 

(2.13) 
for an explicit 



Here {ip'^{t))R stands for the renormalized composite operator [see eq. 
expression] . 

In the N = oo limit the field xi^^t) decouples and does not contribute to the equations 
of motion of either the expectation value or the transverse fluctuation modes. 
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III. UNBROKEN SYMMETRY 



A. Evolution equations in the large N limit 



In this case M|. = iM^jp, and in terms of the dimensionless variables introduced above 



the renormahzed equations of motion are found to be (see references P0|,|2 



fj + 'r] + -rf + g r]{T) S(r) = 



<^,(r) = 



1 



v{0) = Vo , ^7(0) = 



Hence, 



M\t) = 1 + r/(r)2 + g E(r) 



(3.1) 
(3.2) 

(3.3) 
(3.4) 



plays the role of a (time dependent) effective mass squared. 

As mentioned above, the choice of Qq determines the initial state. We will choose these 
such that at t = the quantum fluctuations are in the ground state of the oscillators at the 
initial time. Recalling that by definition gTj{0) = 0, we choose the dimensionless frequencies 
to be 



The Wronskian of two solutions of ( |3.2| ) is given by 

W [<fq, <fq] = 2t , 

while 5'S(t) is given by the self-consistent condition p0| , p7 

roo ( 1 

5S(r) = g q'^dqU ^g(r) ^ - — 



(3.5) 



(3.6) 



Oiq 



2g3 



-v'o+v\r)+g^r) 



(3.7) 



We thus see that the effective mass at time r contains all g-modes and the zero mode at the 
same time r. The evolution equations are then nonlinear but local in time in the infinite N 
limit. 



B. Particle Number 

Although the notion of particle number is ambiguous in a time dependent non-equilibrium 
situation, a suitable definition can be given with respect to some particular pointer state. 
We consider two particular definitions that are physically motivated and relevant as we will 
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see later. The first corresponds to defining particles with respect to the initial Fock vacuum 
state, while the second corresponds to defining particles with respect to the instantaneous 
adiabatic vacuum state. 

In the former case we write the spatial Fourier transform of the fluctuating field ip{x, r) 
in (|2.8|) and its canonical momentum 11 (x, r) as 



Mr) 



V2 



with the time independent creation and annihilation operators, such that annihilates the 
initial Fock vacuum state. Using the initial conditions on the mode functions, the Heisenberg 
field operators are written as 



^,{r)=U-\T) ^,(0) U{t) 



1 



2n 



n,(r) 



U-\r) n,(0)W(r) 




S,(r)=W-^(r) a,W(r) 



with W(r) the time evolution operator with the boundary condition W(0) = 1. The Heisen- 
berg operators dg^r) , aj(r) are related to a^, by a Bogoliubov (canonical) transformation 
(see reference for details). 

The particle number with respect to the initial Fock vacuum state is defined in term of 
the dimensionless variables introduced above as 



(aJ(r)Sg(r)) 



(3.8) 



We consider here zero initial temperature so the occupation number vanishes at r = 0. 

In order to define the particle number with respect to the adiabatic vacuum state we 
note that the mode equations ( |3.2| , |5.2|) are those of harmonic oscillators with time dependent 
squared frequencies 



u''JT) = q^ + l + ri\T)+gEiT). 



(3.9) 



When the frequencies are real (as is the case for unbroken symmetry), the adiabatic modes 
can be introduced in the following manner: 



Mr) 



n,(r) 



20;, (r) 



(3.10) 
(3.11) 
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where now a;g(r) is a canonical operator that annihilates the adiabatic vacuum state, and is 
related to Og, a J by a Bogoliubov transformation. This expansion diagonalizes the instanta- 



neous Hamiltonian in terms of the canonical operators ^(t) 
number is given by 



aHr) 



The adiabatic particle 



(3.12) 



These adiabatic modes and the corresponding adiabatic particle number have been used 
previously within the non-equilibrium context [T^ and will be very useful in the analysis of 
the energy below. Both definitions coincide at r = because u!q{0) = Qq. (For non-zero 
initial temperature see refs. [^,^,0). 

It is the adiabatic definition (|3.12| ) of particle number that will be used in what follows. 



The total number of produced particles N"''^{t) per volume |Mij|^ is given by: 



ad. 



(27r)3 t ^ ' 



(3.13) 



The asymptotic behaviour of the mode functions ensures that this integral converges |20,27 



C. The early time evolution: parametric resonance 

Let us briefly review the dynamics in the weak coupling regime and for times small 
enough so that the quantum fluctuations, i.e. 5'S(r), are not large compared to the 'tree 
level' quantities p0| , |27[| . 



Since S(0) = 0, the back-reaction term 5'S(r) is expected to be small for small g during 
an interval say < r < ri. This time ri, to be determined below, will be called the 
nonlinear time and it determines the time scale when the backreaction effects and therefore 
the quantum fluctuations and non-linearities become important. 

During the interval of time in which the back- reaction term 5'S(r) can be neglected 



eq. (|3.1|) reduces to the classical equation of motion (in dimensionless variables) 

fi + rj + 7f = Q. (3.14) 



The solution of this equation with the initial conditions ( p.3[ ) can be written in terms of 
elliptic functions with the result: 



ri{T) = r]o cn (^T^Jl + r]^, 



(3.15) 



where cn stands for the Jacobi cosine. Notice that rj^r) has period = AK{k)l\Jl + ?7o, 
where K{k) is the complete elliptic integral of first kind. In addition we note that since 

r]{T + 2uj) = -t]{t) , (3.16) 
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1 + ff^T) is periodic with period 2uj. 

Inserting this form for ri{T) in eq. (|3.2| ) and neglecting 5'S(r) yields 



dr2 



+ g" + l + r7o"cn" rVl + r/^,A; 



V9,(t) = 



(3.17) 



This is the Lame equation for a particular value of the coefficients that make it solvable 
in terms of Jacobi functions p0| , p7| . 

Since the coefficients of eq.( |3.17|) are periodic with period 2uj, the mode functions can 
be chosen to be quasi-periodic (Floquet type) with quasi-period 2uj. 



f/,(r + 2a;) = e*^(<?)f/,(r), 



(3.18) 



where the Floquet indices F{q) are independent of r. In the allowed zones, F{q) is real and 
the functions Uq{T) are bounded with a constant maximum amplitude. In the forbidden 
zones F{q) has a non-zero imaginary part and the amplitude of the solutions either grows or 
decreases exponentially. The mode functions v'g(T) obey the boundary conditions eq. (p. 21) 
and they are not Floquet solutions. However, they can be expressed as linear combinations 
of Floquet solutions as follows. 



2ja 



q L 



1 - 



2iVtq 



2iQg 



Uqir) 



(3.19) 



where Uq{0) = 1 and 



= -2q 



\ 



f + 1 + g2 



(3.20) 



We find two allowed bands and two forbidden bands [pOjpTl for eq. ( p.lTp . In the physical 
region > the allowed band corresponds to 



and the forbidden band to 



1^ < g2 < +00 



(3.21) 



(3.22) 



The modes in the forbidden band, < g < r^o/ V^, grow exponentially with time (parametric 
resonance) while those in the allowed band, rio/\/2 < q < oo, oscillate in time with constant 
amplitude. Analytic expressions for all modes were given in p0| , p7| . The modes from the 



forbidden band < g < r^o/v^ dominate S(r). For < r < ri, S(r) oscillates with an 
exponentially growing amplitude. This amplitude (envelope) Se„t,(T) can be represented to 
a very good approximation by the formula [pO| ,^ 



-.Bt 



(3.23) 
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where B and N are functions of rjo given by 

B{7]o) = 8 Jl + vi q (1 - 4g) + Oif) 



4 r (4 + 3 772) /4 ^52 



and the elhptic nome q can be written as a function of r/o as 

with an error smaller than ~ 10"^. 

Using this estimate for the quantum fluctuations S(r), we can now estimate the value of 
the nonlinear time scale ri at which the back-reaction becomes comparable to the classical 
terms in the differential equations. Such a time is defined by g^{Ti) ~ (1 + ?7q/2). From the 
results presented above, we find 



1 

n ~ log 



iV(r^o) (1 + vim 



(3.26) 



The time interval from r = to r ~ ti is when most of the particle production takes 
place. After r ~ ri the quantum fluctuation become large enough to begin shutting-off the 
growth of the modes and particle production slows down dramatically. This dynamical time 
scale separates two distinct types of dynamics, for r < ri the evolution of the quantum 
modes v^g(r) is essentially linear, the backreaction effects are small and particle production 
proceed via parametric amplification. Recall that the zero mode ?7(r), obeys the nonlinear 
evolution equation (|3.14|) . For r > ri the quantum backreaction effects are as important as 



the tree level term ^^(t)^ and the dynamics is fully non-linear. 

We plot in fig. 1 ri as a function of the initial amplitude rjo for different values of the 
coupling g. 

The growth of the unstable modes in the forbidden band shows that particles are created 
copiously (~ 1/(7 for r ~ ri). Initially, (r = 0) all the energy is in the classical zero 
mode (expectation value). Part of this energy is rapidly transformed into particles through 
parametric resonance during the interval < r < ti. At the same time, the amplitude of the 
expectation value decreases as is clearly displayed in fig. 2. We plot in fig. 3 the adiabatic 
number of particles [as defined by eq.( p.l3| )] as a function of time. 

The momentum distribution of the produced particles follows the Floquet index and is 
peaked at g ~ |?7o {1 — q) [^,0] this is shown in fig. 4. 



IV. ASYMPTOTIC NONLINEAR EVOLUTION 

A. Numerical Analysis 

In the previous section we have summarized the dynamical evolution in the linear 
regime in which the backreaction effects can be neglected and estimated the first new, non- 
perturbative dynamical time scale ri as that beyond which the dynamics is fully non-linear. 
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In this section we present the time evolution after the nonhnear time ri. That is, when 



the backreaction 5'S(r) is important and the full solutions to the non-linear equations ([O 
|3.3| ) is needed. We have implemented a refined numerical treatment for a wide range of 
initial amplitudes and couplings. The numerical method uses a fourth order Runge-Kutta 
algorithm and 16-point Gauss integrations for the integrals over q and is appended with 
a FFT analysis to determine the frequency spectrum of the oscillatory component. The 
precision of our results is better than one part in 10^. 

To begin with, we observe that 5'S(r) and //^(r) oscillate with the same frequency and 
opposite phase. Thus, a remarkably cancellation takes place between these two terms in 
the effective mass squared. This phase opposition is analogous to Landau damping |]T3[. 
One sees such cancellation comparing fig. 2 for ^^(t), fig. 5 for 5'S(r) and figs. 6 for A1^(r). 

Moreover, we see that A1^(r) tends to a constant value for r oo. We find numerically 
that this value turns out to be 

-ML = 1 + 1 (4.1) 

for the values of g and rjQ considered in figs. 1-12. (up to corrections of order g that are 
beyond our numerical precision). It must be noticed that coincides with the lower 

border of the allowed band ( |3.21| ). 

Furthermore A^^(r) approaches its asymptotic limit (|4.1|) oscillating with decreasing 
amplitude. More precisely, using a detailed numerical analysis of the asymptotic behavior 
and fast Fourier transforms we find from our numerical results for r > ri [see figs. 6a and 
6b], 

M\r)=Ml^ + ^-^ + o[^^) (4.2) 

with 

pi(r) = Ki cos[2 Moo T + 2a2 log(r/ri) + 71] + K2 cos[2 MqT + 262 log(r/ri) + 72] , (4.3) 

where i^i, i^25 7i and 72 are constants and an excellent numerical fit for the coefficients a2 
and 62 is given by 

02 ^ 0.16 In - + 0.6 

9 

62 0.6 - 0.16 In- (4.4) 
9 

within a wide range of (weak) couplings and initial values of //(O). We also find that the 
coefficients Ki, K2 vary linearly with In(^) a result that will be obtained self-consistently 
below. 



1. The evolution of the expectation value and mode functions 

Since the effective mass tends asymptotically to the constant value A4oo, the expectation 
value rj^r) oscillates with frequency A4oo and the g-modes (Pqir) with frequency 
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oo{q)^^q' + Ml, . (4.5) 

[Notice that uj{q) = ujg^r = +00)]. These oscillation frequencies are confirmed by the 
numerical analysis of the evolution of the expectation value and the g-modes. Figs. (7 
a-c) display the momentum distribution of the created particles at different times. One of 
the noteworthy features is that whereas up to time t ^ ti the distribution only has one 
peak at the value of maximum Floquet exponent, for larger times the backreaction effects 
introduce new structure and oscillations, keeping the borders of the band fixed throughout 
the evolution in the non-linear regime. 

The numerical results displayed in figs. (7 a-c) show that the position of the main peak 
go(T), decreases with time. We performed a numerical fit for the time dependence of the 
peak position and found its behavior to be well described by the estimate: 

qlir) ^ ^ . (4.6) 



with the constant Ki introduced in equation (|4.3| ) above. We will provide an analytic, 
self-consistent description of this behavior below. 

g-modes above and below go behave quite differently. Modes with q > qo{T) oscillate 
in time with constant amplitude. Modes with q < qoir) also oscillate but with increasing 
amplitude. 

Eq. ( |4.6| ) shows that the peak position in q decreases monotonically as ~ 1 / ^/T. As time 
evolves, more and more g-modes cross the peak and become purely oscillatory. Only the 
amplitude of the g = mode (which is not to be confused with the expectation value rj^r)), 
keeps growing. As we shall discuss in detail below, there is a band of non-linear unstability 
for < g < qoir). We also find numerically, that a second non-linear resonance band appears 
just below g = -/^o/v^ for gi(r) < g < rjo/\/2, and we find numerically that 

« I - ^ . (4.7) 

Z T 

However, the growing modes in this upper band give a much less important contribution to 
the physical magnitudes than the first band. 

The modes in between, qoi^r) < q < qiij) , oscillate for times r > ri with stationary 
amplitude Mq{r). 

This crossover behaviour of the modes can be expressed by introducing a g-dependent 
time scale beyond which the modes become oscillatory. Such scale is given by 

rM) = ^ (4.8) 

for the lower nonlinear band and 

miq) = (4.9) 
f -g2 

for the upper nonlinear band. 

Both ipq=o{T) and 77 (r) obey the equation (|3.1| )and are linearly independent solutions, 
their difference arising from the initial conditions. ipq=o{T) has growing amplitude while the 
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amplitude of ?7(r) decreases with time. Since these are hnearly independent solutions of the 
same equation their Wronskian is a non- vanishing constant. Therefore if one solution grows, 
the other independent solution must decrease in order to respect the wronskian condition 
( p.6|) . Since the total energy is conserved rj^r) must necessarily be a decreasing solution. 

The fact that the zero- mode amplitude 7]{t) vanishes for r = oo implies that all the 
available energy transforms into particles for r = oo. This conclusion which will be fur- 
ther clarified in what follows is a consequence of the non-linear dynamics. It is the more 
remarkable because the particles produced are massive and therefore there is a threshold 
to perturbative particle production. It will be seen in detail below that the particle pro- 
duction in this regime is a truly non-perturbative phenomenon associated with non-linear 



resonances. 



For r > Ti,Tj{q),Tii{q) the effective mass squared tends to a constant [see eq. ( |4.1|) ], 
therefore, the asymptotic behavior of '^qir) is given by 

^g(r) "=°" Ag e'"(«)" + Bg e-^'^")" + O (^^) (4.10) 

with uj{q) given by eq.( [4.5|) . It is then convenient to define the following functions 



-iLu(q) T 



uj{q) 



<^g(r) + 



uj{q) 



which for t > Ti are slowly varying functions of r, with the asymptotic limits 



lim Aq{T) = Ag ; lim 5q(r) = Bg 



(4.11) 



(4.12) 



We can thus express the mode functions (Pg^r) in terms of Agi^r) and -Bg(r) as follows, 

(^,(r) = Agir) e'"(«)^ + 5,(r) e-'"(^)" . (4.13) 
We obtain from eq.( [4.11| ) for the square modulus of the modes, 

|^,(r)|2 = \AgiT)\' + \BgiT)\' + 2\AgiT) BgiT)\ COS [2^(g) r + </.,(r)] (4.14) 
where we have set 

(4.15) 



Mr) Bg{T) 



\Aq{T) 5,(r)|e^^«W 



The wronskian relation ( |3.6| ) implies that the functions Ag{T) and Bg^r) are related asymp- 
totically through 



mi 



\AJt) 



uj{q) 



(4.16) 



plus terms that vanish asymptotically. The virtue of introducing the amplitudes 
Ag{T) ; Bg{T) is that their variation in r is slow, because the rapid variation of the mode 
functions is accounted for by the phase. 
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Figs. 8-11 show the (scaled) modulus , 



(4.17) 



and </>q(T) for some relevant cases. As shown in these figures Mq(r) and (f)q{T) do not exhibit 
rapid oscillations with period 27i/uj{q) and 27r/A^oo which are present in (pg{T) and //(t), 
respectively. That is as anticipated above, Mg(r) and (pqir) vary slowly with r. 

For small coupling g, \(pq{T)\'^, |i?q(r)p and |Ag(r)p are of order 1/(7 for g in the forbidden 
band and times later than ti ||20|j27[] . Therefore, Mq{T) becomes of order one after the non- 
linear time scale for modes inside the band, and is perturbatively small for modes outside 
the band. Moreover, eq. (|4.16|) implies that |i?g(r)p = |y4g(T)|^[l + 0(g)] and for modes 
inside the band we can approximate eq. ([4.14| ) as follows. 



Mqir) 



{1+ cos [2^(g)r + 0,(r)]} [1 + 0(^7)], 



(4.18) 



for < g < r]o/V2. This expression is very illuminating because it displays a separation 
between the short time scales in the argument of the cosine, and the long time scales in the 
modulus Mg(r) and phase (pqir). 

We now introduce slowly varying coefficients for the order parameter ^^(t). Let us define 



Dir) 



2 

00 



-M.00 T — arctan 



fjir) 



(4.19) 
(4.20) 



Using the result that asymptotically the effective time dependent mass reaches the asymp- 
totic limit A4oo and the fact that 77 (r) is a real function we write 



rj{T) = D{t) cos [Moo t + 0(r)] 

We plot in figs. 12 D[t) and 0(r) as functions of r 



i + o[- 



(4.21) 



We begin our numerical analysis by considering the g = mode function. After an 
exhaustive analysis we have found that Aq^oij) and -Bg=o(T) exhibit power behaviour for 
r > Ti [see figs. 8-10]. To our numerical precision these power laws can be fit by the 
following form 



A=o(r) 



Bq=o{r) 



where the numerical results yield for the anomalous dynamical exponents 

ai ^ 0.27 



(4.22) 



(4.23) 



while a2 is the same as in equations 

The behaviour ( [4.22|) appears also in the evolution of the expectation value ?7(r) but 
with the growing power of r absent [ Ci = C[ = ) resulting in that the order parameter 
decreases with time exhibiting a logarithmic phase: 
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r]{T) = Do [ — ) cos[A^oor + a2log(r/ri) + /o] 



l + O 



(4.24) 



where /o is a small constant. Therefore, comparing with eg. ( [4. 21 ) we find the remarkable 
result that the amplitude of the expectation value relaxes with a dynamical power law 
exponent: 



D(r) = Dn 



and a logarithmically varying phase 











1 - 


-(7)] 









(r) = a2log(;^)+/o + o(^ 



(4.25) 



(4.26) 



The time Ti appears here since it is the natural time scale for the non-linear phenomena. 

The q 7^ 0-modes also grow with time with a power-like behaviour for < g < qoir) 
but with a larger power than the q = mode. [See figs. 8-10]. Such growth is definitely 
milder than the exponential increase of the modes inside the forbidden band in parametric 
resonance. Our interpretation of this phenomenon is that < q < qolr) is a non-linear 
resonant band. It is not a resonance in a linear differential equation (as it is parametric 
resonance), but a new nonlinear effect which is a consequence of the backreaction of the 
quantum fluctuations through gT,{T). 

A second non-linear resonance band appears just below q = rio/ \/2 for qiir) < q < ^J2. 
However, we flnd numerically that the contribution from this upper band to the physical 
quantities such as particle production, is much smaller than the flrst band. The modes in 
between, q^ij) < q < qii^r) , oscillate for times r > ri with stationary amplitude Mg(r). 

The nonlinear resonant bands become narrower as a function of time, i.e. qo{T) and 
r7o/\/2 — ?i(t) decrease with time (qiir) increases). 

The growth of the amplitudes Mq{T) in the nonlinear resonant bands for a flxed q stops 
when q crosses the borders qoir) or qiir). After that time, such g-modes oscillate with 
constant amplitude this behavior is displayed in flg. 8d. There is a crossover for q ~ qo{T) 
and for q ~ qi{'r) from monotonic growth to oscillatory behavior. 

The phase ^^(t) exhibits an analogous behaviour (see flg. 9d). 

The particle distribution exhibits marked peaks at g ~ Qo{t) and at g ~ gi(T), which 
are clearly displayed in flg. 7. Notice that the peak near g qiif) has a much smaller 
amplitude. 

The mode exactly at g = r/o/v^ has an analogous behaviour to the g = mode. [Compare 
flgs. 13 with flgs. 8]. We flnd to our numerical accuracy that the amplitudes behave as 



^q=ryo/v^(^) 



r 



ib2 



El t'' + E2 T-'^ 
The numerical calculations yield for the dynamical exponents 

61 ^ 0.19 

and 62 is the same exponent as in eqs. (^4.3|) - (^4.4|) . 



E[t^' + E'^ T-" . (4.27) 



(4.28) 
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The growth of the g-modes in both nonhnear resonant bands leads to particle production. 
We see from fig. 3 that the number of particles continues to grow after the nonlinear time 
r = Ti. Although this growth is much slower than before t = ti, the total number of 
particles produced after the time ti is substantial and turns to be of the same order of 
magnitude than those produced before ri. 

The adiabatic number of produced particles N^'^{t), can be expressed for late times in 
terms of the mode amplitudes Mq{T) as follows: 

Nfir) = ^ uj{q) M,(oo)^ + O {^^ + O (g) (4.29) 



where we used eqs.( p.l2D , ( |4.1CI| ) and ( |4.17| ). Notice that asymptotically the adiabatic particle 



r Vr^ 



number depends solely on the long time scale as the terms containing the fast oscillating 
function cos [2 uj{q) t + (pqir)] cancels out to order for large r. This is one of the important 
advantages of this definition of the particle number. 

For T >> Ti the total number of produced particles approaches its asymptotic value 
A/'"'^(cx)) as 

Af'^r) = Af^^'ioo) -^ + (^) , (4.30) 

where 

1 fVo/V2 

?A/""(oo) = - 
and G is positive. 

The numerical analysis show that gAf°''^{oo) and gG depend very little on g for small 
g < 10~^. Both J\f°''^{oo) and G grow with r]Q. Precise numerical fits yield the behavior 

gAf'^oo) ^ 0.007 r]^-^ (4.32) 

for a wide range of couplings and rjo. 

At this point we summarize the results from the numerical analysis for the unbroken 
symmetry case: 

• The effective time dependent mass reaches a finite asymptotic value A^oo in the form 



gX^'^(oo) = —J^ dq uj{q) Mqioof + O (g) (4.31) 



given by equations ( |4.2| ,0)-(0). This in turn means that the modes become free 
asymptotically with plane wave behavior and the non-linear self-consistent coupling 
between modes vanishes. The large N limit yields free modes in the infinite time limit. 

For weak coupling and for r > ri there is a separation of time scales, with a short 
time scale corresponding to the oscillations with frequencies corresponding to a mass 
~ A^oo and a longer time scale that depends on ti. 

The amplitude of the expectation value relaxes with a power law with non-universal dy- 
namical exponents and logarithmic phases that vary solely on the long time scale. The 
expectation value vanishes asymptotically, despite the fact that its energy is dissipated 
into massive particles for which there are perturbative thresholds for production. The 
relaxation mechanism is non-linear and clearly non-perturbative even at long times. 
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• For r > Ti there are non-linear resonant bands which form at the edges of the origi- 
nal band for parametric amplification. The width of these non-linear resonant bands 
vanishes asymptotically resulting in that all modes oscillate harmonically for asymp- 
totically large time. For r = oo both unstable non-linear bands { < q < QqIt) 
and Qiir) < q < r^o/v^ shrink to zero). The crossover from power-like to oscillatory 
behaviour takes place at the g-dependent time scales given by eqs.( [4.8|) -( [4.9|) . 

• The particle distribution N^'^{t) has a finite and nontrivial limit for r oo. In 
particular, a consequence of the non-linear resonant bands is that N^'^{oo) will be 
peaked at g = 0. The asymptotic form of the distribution is a function of the initial 
conditions and the coupling g. In particular, N^^{oo) is of order 1/g for q < r/o/v^ 
and it is of order one for q > r^o/v^- That is, the support of the particle distribution 
valid for short times r < ri survives for all times including t = oo. Furthermore, 
for weak coupling the large number of particles inside this band allows us to interpret 
this asymptotic state as a non-perturbative semiclassical condensate in the unbroken 
symmetry phase that has formed dynamically through the relaxation of the initial 
energy. 



B. Asymptotic Analysis I: Perturbation Theory 



In the previous section we presented an exhaustive numerical study of the evolution 
of the mode functions and the expectation value. In this section we provide an analytic 
perturbative approach to explain and understand the numerical results. 

In order to study analytically the asymptotic behaviour for late times, it is convenient 
to write the equations for the expectation value and mode functions as follows: 



rf2 

— + Ml + wir) 
(f 

— + MI, + w{t) 



^,{r) = 
r/(r) = 



where 



w[t] 



r 



M 



O 



(4.33) 



(4.34) 



and Pi(t) given by eq. ( ^^ ) will be treated as a small perturbation for r >> ri. 

These equations can be written as integral equations using the proper Green's function. 
That is. 



, sinci;(g)(r' — r) 



dr' 



uj{q) 



(4.35) 



Here we used the advanced Green's function that obeys 
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dr 



[T - r 



sin[a;(g)(r' — r) 



uj{q) 



(5(r - r') . 



(4.36) 



Since w{t) = O (^^^ we can generate the asymptotic expansion for V5q(r) just by iterating 
eq.(|;35D. We find, 



uj{q) 



+ o 



(4.37) 



The integrals here can be performed in closed form up to terms of O (^^^ by using eg. ( |4.3| 
for Pi{t). The result is given by 



A. 



i^i sin\I^i(r) ^ sin \I^2(i" 



8uj{q) T 
+ K2 



AiM^uj{q)T AiMoUj{q)T 

(r)-2 T(^(q)) -i(<^i(r)+2 r a;(g)) 



uj{q) -Mc 



uj{q)+Mc 



gj(<I'2(T)-2ru;{g)) g-i(1'2(T)+2 r 



+ 



+ [Bg 

A. 

+ 



uj{q)~Mo u{q)+Mo 
Kisin\I^i(r) i^2sin\I^2(T 



4:iMooOo{q)T AiMoUj{q)T 



\uj{q) T 



g-iC-I-l {t)-2 T Lu{q)) gj(<I'i (r)+2 r Lu{q)) 
+ 



Lj{q)-Mc 



uj{q)+Mc 



+ K2 



+ o{- 



-j(1'2(r)-2r<^(g)) 

u{q) - Mo 



+ 



,i(*2(T)+2T<^(g)) 

uj{q)+Mo 



(4.38) 



where 



and 



^i(r) = 2A^oor + 2a2log- + 7i 



^2(7-) = 2;V/lo ^ + 2&2log — + 72 • 



These expressions display resonant denominators for u^q) = Moo and u^q) = Mq. These 
resonances correspond to g = and q = r]o/\/2, respectively. This perturbative approach 
is expected to be valid when the first order correction is smaller than the zeroth order. A 
necessary condition for its validity is given by 



uj{q) T [io{q) - M 
This implies for q significatively smaller than M 



< 1 . 



(4.39) 
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(4.40) 



where, we approximated 



Thus in the regime where eq. ( [4.40|) holds the behavior of the mode functions is oscillatory 



and given by (|4.38|) . This is in agreement with the numerical results discussed in sec. IV. A. 
The equality sign in eq. ( |4.40| ) yields the eq. ( f4.6|) for the peak position which was found 



numerically and therefore now interpreted as the result of a resonance condition. Such 
peaks can be seen also in the mode functions amplitudes displayed in figs. (8,11). 

The resonance at g = ?7o/a/2 can be treated analogously. The necessary condition for 
the validity of the perturbative approach is then 



uj{q) T [uj{q) - Mc 
Near q = r^o/v^ we can write, 

V2Mo 

Therefore, the oscillatory behaviour ( [4.38[) applies for 



^-q>^. (4.42) 



The equality sign in eq. ( |4.42| ) yields the peak positions of the mode functions amplitudes 
near such resonance: 

gi r = ^ . 4.43) 

These results are in remarkable agreement with our numerical calculations [see figs. 8,11]. 

The position of the main peak in the particle distributions precisely corresponds to the 
situation where q^ is balanced by the amplitude in the 'potential' Namely, for 

q^ > Ki/t, we have oscillating modes and for < Ki/t, resonant (growing) modes. The 
same argument applies to the secondary peak. 



C. Asymptotic Analysis II: Multitime Scales 

The perturbative analysis took us a long ways towards understanding the presence of 
the non-linear resonances associated with the backreaction effects and revealed the position 
of these resonances in complete agreement with the numerical study. 

However to describe the evolution of the modes inside these bands the perturbative 
approach is insufficient and a non-perturbative method of resumming the potential secular 
terms associated with the resonances must be implemented. 
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The main observation from the numerical analysis is that for weak coupling there are two 
widely separated time scales, the short time scale associated with oscillations with frequency 
determined by the asymptotic value of the effective mass, and a long scale associated with 
the non- linear time ti. This suggests to implement a multitime scale analysis p8[ which 
resums the secular terms and results in a uniform expansion. This method implements a 
dynamical renormalization group resummation which was already implemented successfully 
in non-equilibrium evolution in quantum field theory and previously applied to quantum 
mechanical problems 

In this section we implement the method of multitime scales to the equations for the 
expectation value and the q — modes which we write in the following form 

dr^ 



+ q^ + M\T) 



V9,(r) = 



(4.44) 



r]{r) = 



(4.45) 



and use the asymptotic behavior of the effective mass squared obtained from the detailed 
numerical analysis A1^(r) for r > ri 



Ml 



= Ki cos [2Moor + 2a2 ln(r/ri) + 71] 
+ K2 cos [2MoT + 262 ln(r/ri) + 72] 



(4.46) 
(4.47) 



with Ti ^ In i [see eq. ( |3.26| )] being the nonlinear time scale. 

As emphasized above, the non-perturbative dynamics has generated a new time scale ti 
and for weak coupling there are at least two widely separated time scales, the short time 
scale corresponding to oscillatory behavior of 0{Ai^) and the long time scale for non-linear 
relaxation of the order of ri. 

In order to implement the multitime scale analysis it is convenient to introduce the small 
quantitity e and the following two time variables (Tq and Ti) by 

e = — ; To = r ; Ti = eTo = — 
n Ti 

-^ = Do + eDi ; = ; n = 0, 1 

ar din 

and to write Pi{t)/t in a manner that displays at once the dependence on the short and 
long time scales 

= e r(ro,Ti) 

T 

r(To, Ti) = ^ COS [2M00T0 + 2a2 In(Ti) + 71] 
1 1 

+ ^ cos [2M0T0 + 262 In(Ti) + 72] (4.48) 

-L 1 
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To C(e) the multitime scale analysis of the asymptotic time dependence of 77 ; ipq begins 
by proposing the following uniform perturbative expansion for the solution 

r/(To, Ti) = r/W(To, Ti) + e r7«(To, T,) + ■ ■ ■ 
^g{To, Ti) = <)(To, Ti) + e yp«(To, Ti) + ■ ■ • (4.49) 



D. g = modes: rj and 93^=0 

We generically call /(To,Ti) both 77 and iPg=o. The only difference between these is that 
whereas rj is always real, (pq=o is complex, this difference will be accounted for in the final 
form below. Comparing powers of e, we find the following equations for the g = modes to 
first order in e 







[2DoDi + r(ro,ri)]/(°)(ro,T, 



The solution to (^4.501) is obviously 



(4.50) 



(4.51) 



where for rj the reality condition implies B{Ti) = A*{Ti). If the solution of ( [4.50 ) is sought 
in terms of the Green's function of the differential operator on the left hand side, one finds 
that the term proportional to cos [2 00 ^o] in T(Tq,Ti) would give rise to secular terms. 
Therefore the condition for a uniform expansion requires that the coefficients of these secular 
terms vanish. This leads to the following differential equations for the dependence of the 
coefficients on the long time scale Ti, 



DiA- 
We find the solutions, 



gi2a2 ln{Ti)+i7i 
g-i2a2ln(Ti)+i7i ^ _ Q 



(4.52) 



A(Ti) = a± e5[2'^2ln{Ti)+7l] J^iai ^ 5( J^^) = 5^ g-f [2a2 lnm)+7l] J^±ai 



ai 



tan 5 



( ^1 



(4.53) 



ai 
a2 



This solution confirms the power law relaxation found numerically and provides the 
consistency condition 
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This condition is verified numerically to our level of precision. In addition for weak 
coupling the numerical evidence gives << (see equations (|4.4| ) and (K:.23|) leading to 



(4.54) 



The final form of the solution is given by a linear combination of the two independent 
solutions above, yielding 



/(ro,ri)=Tf c+ cos 

+ T^"^ C. COS 



7WooTo + a2ln(ri) + ^ + ^ 



A<ooTo + a2ln(Ti) + 



7i 



0{t 



where for 77 the coefficient C_ is real and C+ = 0, whereas they are complex and both 7^ 
for (pq=o. The 0{e) correction quoted above is hound as a function of time as it arises from 
the perturbative solution without secular terms. 

Therefore we quote the final form of the solutions 

-11 r s 

cos MocTq + a2 In(Ti) + y - 2 



r/(r) = 1^0- 



ai 



COS 



M^n + a2 In(Ti) + y + ^ 



COS MooTq + 02 In(Ti) + y - 2 



In addition, we have checked the constancy of the Wronskian 



C.C* - CIC. 



leading to the conclusion that neither of the coefficients C± can vanish for ipq=o- It is a 
matter of straightforward algebra to find that these are indeed solutions of the equations 
( [4.44| ) (for q = 0) and ( [4.45| ) up to terms that fall off faster in time. Furthermore, a uniform 
perturbative expansion in e can now be carried out to the next order. 



E. g 7^ modes 



Proposing a uniform e- expansion for the mode functions and keeping only up to 0(e) 
we find the equation 



D'i + 2eDoD, + + Ml + eT{To, T,) (^("HTo, T^) + e ipf\T,, T, 







(4.55) 



Since F contains the oscillating factors cos [2A^oo^o] aiid cos [2A^o^o] a naive perturbative 
expansion in e will produce secular terms for ^ A4oo ,-M.o, i-e- for g ^ , A4l — 
Ail = ?7o/2, these were the values for which the resonant denominators in the perturbative 
expansion vanish. Therefore for these values of q we must implement a multi-time scale 
analysis to resum the secular terms. 
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1. Small hut nonzero q 

A consistent expansion in e can be implemented by writing q"^ = e qf with qi ^ C>{1). 
This is a nonrelativistic approximation since then is much smaller than A1^. The zeroth 
order solution is clearly 



Secular terms in a naive perturbative expansion will arise from the term proportional to 
Ki in Pi{t). 

It is convenient to define the coefficients 



(4.56) 



Requesting that the coefficients of the secular terms in the perturbative solution vanish we 
obtain the following differential equations 



Dibg- 



ia2aq 
ia2hq 



iql iKi 

a„ - ... „ ba = 



iql 



a„ = 



These equations simplify considerably by introducing the variable z as 

_ q^T, 



2Mc 



(4.57) 



Then, eqs.( |4.56| ) can be rewritten as two decoupled second order differential equations: 

72 J „2- 



d d . ^ a\ 

z—^ + — + z - i- 2a2 

dz'^ dz z 

d'^ d a\ 

Z-— + — + z + i ~2a2 

dz"^ dz z 



ctqiz) = 



bg{z) = 



(4.58) 



These are conffuent hypergeometric equations with the solutions: 



Mz) = z^'^--^ M,,^^,^{2tz) , B, 



z = z 



5 M 



ia2- 



(4.59) 



where Mx,^{z) stands for a Whittaker function [^. We find the asymptotic behavior to be 
given by |^ 



B,{z) (2f 



and 



Ag{z) 
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.1 r(l±2ai) 
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{2^] 



ia2- 



1 r(l±2ai) 

r(l + ia2 ± di] 



In figs. 15 we plot and the phase of e~*^ Aq{z) as a function of z. One sees that 

the behaviour of the numerically computed modes in figs. 8c, 8d, 9c and 9d is accurately 
reproduced. 



Using the integral representation for the solutions [Q, we can find where the functions 
Aq{z) and Bq{z) oscillate with z and where they do not. 

Aq{z) = k z''''^'^' f - e42-*+'^2i°g[¥l} [t(i _ t)f-^ 

J t 

where A; is a constant. This integral has stationary points at: 

1 



t 



1 ± Wl - 



2a2 



For z > 2a2 the stationary points are real indicating an oscillatory behaviour whereas 
they are complex for z < 2a2 implying a non-oscillatory behaviour. We see from eqs.( [4.57D 
and (|4.54| ) that the condition z = 2a2 precisely corresponds to the peak position ( [4.61) . 



2. q^^Ml-Ml 

In this region of momentum we write = Ai^ — M.^ + eg| in order to implement a 
multi-time scale analysis. In this case, secular terms in the perturbative expansion in e will 
arise from the term proportional to K2 in Pi(t), i.e. the term cos[2A^o^o]- 

Using the same notation as in the previous subsection, the zeroth order solution is now 

(^f (To, Ti) = ^,(Ti) e*^°^° + i3,(Ti) e"*^"^" (4.60) 
Defining the coefficients 

a^(Ti) = e-^l2^2WT0+72]^^(j.^) 

Requesting that the coefficients of the secular terms in the perturbative solution vanish we 
obtain the following differential equations 



Dittq + ih2 

Dibq - ib2 



aq 


■ ql 


. K2 


Ti 








K2 







As in the previous case these coupled first order differential equations reduce to two 
decoupled confluent hypergeometric equations with solution: 

Aq{y) = y''-^~'^ M,,^_.^^,piy) , S,(y) = y^^-^ M_,,^_ (2zy) 
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where 

_ ('toV2 - q') 

The transition from non-oscillatory to oscillatory behaviour takes place here aX y = 262 
implying that the modes with 

^2 T 

will grow in time. 

The solutions for the mode with q = r]o/^/2 correspond to exponents 



2 



This analysis leads to the conclusion that as a consistency condition, the constants Ki 
and K2 can be expressed as follows, 



i^i ^4 ^^1 + 77^/2 (0.16 log ^-^ + 0.6) , ~ 4 1/1 + r/^ (-0.16 log ^"^ + 0.6) , (4.61) 

where we used eqs.(U), ( ^:23D and dOSj ). 

Therefore the effective mass squared behaves as 

4 

M\r) = Ml^ + - {Moc a2 cos[2 Moo r + 2a2 log(r/ri) + 71] 
r 

+ Mob2 cos[2;WoT + 262log(r/ri)+72]} + C»(^) 
We have confirmed these results numerically within our precision. 

These results are noteworthy, by implementing a multitime scale analysis which is the 
dynamical equivalent of a renormalization group resummation p9[ we have obtained power 



law relaxation for the expectation value with non-universal dynamical anomalous dimen- 
sions. The logarithmic phases are clearly a consequence of the 1/r fall off of the potential 
w{t) in the mode equations ( |4.33|) just as in the Coulomb problem. The power laws orig- 
inate from the resummation of the secular terms arising from the non-linear resonances, a 
non-perturbative result. 

F. Energy and Pressure 

The energy-momentum tensor for this theory in Minkowski spacetime is given by 



(4.62) 



Since we consider translationally as well as rotationally invariant states, the expectation 
value of T^" takes the fluid form 
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E 



NV P(r) = < T^\x) >=< T^^{x) >=< T^^x) >=< ^ (V^)^ + (jf - T°°(x) > 



with all off-diagonal components vanishing. 
Hence, 

1 1 ^ 



Pi 



E 



NV 3 
takes a particularly simple form. 

Both E and -P(t) + E can be expressed in terms of the zero mode and the g-modes. 
All derivations including the renormalization procedure can be found in refs. 0,^,^]. We 
just quote the final results in the unbroken symmetry case referring the reader to the above 
references for details. 
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(4.63) 



where 20,27 



q'dq u;,{T)Nf{T 



l + ril 



1 + log 



-4(1 + %') ^^M- 



Since energy is conserved (as can be verified explicitly by using the equations of motion), 
it is equal to the initial value, given by 
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2|M, 



-Rl 



2|M, 



-Rl 



1 + -^ni 



The renormalized energy plus pressure takes the form pD| , p7| , 



2|M, 



R\ 



[e + p{T)] 



2|M, 



-Rl 
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-R 



Jo 



— Q 
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M\t) e{q-K) cP 
3g ^ 12 g3 dr^ 



(4.64) 



For times after ri we can restrict ourselves to the contribution from the expectation value 
•T] and the modes in the band < g < 770/ v^- Modes with q > r]Q/\/2 only yield perturba- 
tively small corrections 0{g). Using eq. ( [4.11| ) we can write the integrands in eqs.( [4.63"D -( [4.65| ) 
as follows: 

9Mt)\' = M,{rf {1 + cos [2 uj{q) r + 0,(r)]} [1 + 0{g)] , 
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g\^,{T)\^ = uj{qf M,{Tf {1 - cos [2uj{q) r + </),(r)]} [1 + 0{g)] , 
Inserting these expressions in eqs.( |4.63| )-( |4.65| ) yields 

e = lr,' + l{l + v')M'ir)-^^^^^^ + fj'^\' dq [q' + M\oo)\ M,{rf + 0{g) . 



Taking now the t ^ oo hmit yields, 



1 , fno/V2 



q'^dq q'' + M\oo) M,(oo)2 + C(^) . 



(4.65) 



We analogously find for e + p(r) 



e + p{T) 



■no 



/V2 



q^ dq 



MJt? 



■no 



/V2 



q^ dq cos [2 uj{q) r + 0g(r)] 



M,{rf + 0{g) 



For large r the integral containing the oscillating cosinus dies off. We thus obtain combining 
both expressions: 



1 /•r)o/^/2 1 

p(oo) = - / qUqM,{^f + -4 + 0{g). 
O JO lo 



(4.66) 



G. Sum Rules and the Equation of State 



Although we do not know the analytic form of the particle distribution for late times 
[see eq. ([4.29| )], we are able to compute its first two moments in the following way. 

First, we can express the quantum fiuctuations g^ij) in terms of the modes using 
eqs.(|;3) and (|1|), 

^S(r)= / q^dqM,{Tf{l + cos[2uj{q)T + (l,,{T)]} + 0{g). 
Jo 

For large r the integral containing the oscillating cosinus dies off. Using now that ri{oo) = 0, 
A4^(r) = 1 + ?7(r)^ + g S(r) and eq.(|4.lD we obtain the first sum rule 



rvo/v2 1 
/ q'dqM,{^)' = -r^l + Oig). 
Jo 2 



Furthermore, equating the expression for the energy for r = oo [eq.( |4.65| )] with its initial 
value [eq. ( [4.64| )], yields the second sum rule: 



V0/V2 1 

qUqM,{ooy = -r^', + 0{g). 
lo 



Combining the sum rules with the expressions for the energy and pressure, eqs.( [4.65"D 
and (|4.66| ) yields 

p{^) = ^Vo + 0{g). 
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and 

We see that the bath of produced particles does not behave asymptotically either as radiation 
or as nonrelativistic matter but their equation of state interpolates between these two limits 
as a function of the initial amplitude of r^. 

For large r^o, we find as expected, radiation behaviour, 



p(00) ^o^oo 



For small r^o, we find a cold matter behaviour. 

We notice that the energy and the asymptotic pressure can be expressed in a form that 
suggests a two component fluid formed by nonrelativistic and massless particles: 

^ and ^ may be interpreted as the contributions of massless and massive particles to the 
total energy. 

The pressure approaches its asymptotic limit p(cxd) oscillating with decreasing amplitude. 
We find from our numerical analysis [see figs. 14a and 14b] that this amplitude falls off as 
~ 1/r. More precisely we find for times later than ri, 

V{r) = p{oc) + ^ + O (1 
r Vr^ 

Here, g'i(r) oscillates with time with the same frequencies 2Aioo and 2Aio as the effective 
mass squared [see eq.(|4.3|)]. 



V. BROKEN SYMMETRY 



In the case of broken symmetry = — |M^| and the field equations in the N = oo 



limit become 20,27,34 



fj — Tj + rf + g Tjij) S(r) = 



V9,(r) = 



(5.1) 
(5.2) 



where S(r) is given in terms of the mode functions ^q{T) by the same expression of the 
previous case, ( |3.7| ). Here, A^^(r) = — 1+ ri{T)'^+g S(r) plays the role of a (time dependent) 
renormalized effective mass squared. 
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The choice of boundary conditions is more subtle for broken symmetry. The situation of 
interest is when < ?7o << 1, corresponding to the situation where the expectation value 
rolls down the potential hill beginning very close to the origin. The modes with < 1 — r^Q 
are unstable and thus do not represent simple harmonic oscillator quantum states. Therefore 
one must chose a different set of boundary conditions for these modes. Our choice will be 
that corresponding to the ground state of an upright harmonic oscillator. This particular 
initial condition corresponds to a quench type of situation in which the initial state is evolved 
in time in an inverted parabolic potential (for early times t > 0). Thus we shall use the 
following initial conditions for the mode functions: 



^.(0) = 




= \lq^ + ^ + vl for <ql = l-r]l 

^g = - 1 + for (f><il ; < r^o' < 1 • (5-4) 

along with the initial conditions for the expectation value given by eq.( |3.3| ). Furthermore 
because the adiabatic frequencies cannot be defined for the modes in the spinodal band, we 
use the definition eq.(|3.8|) for the particle number. 



A. The early time evolution: spinodal unstabilities 



As in the unbroken case, for << 1 we can neglect 5'S(r) in eqs. ( |5.1| - |5.2| ) until the 
spinodal time ti defined to be the time scale at which the quantum fluctuations become 
comparable to the 'tree level' terms. In addition, when the initial value of rjQ is zero or much 
smaller than one, we can neglect 77(r)^ in the mode equations, which simplify to 







(5.5) 



The solution of eqs.( |5.5| ) for the initial conditions ( |5.3| ) takes the form 
1 



q 



i[l + q' 



+ 



1 - + i(l + 



for < g^ < 1. The modes with higher wavenumber g^ > 1 do not grow and their contribu- 
tion to 5'S(r) is subleading at long times. 

This solution exhibits the typical exponential growth of spinodal unstabilities. Very 
soon, the decreasing exponential can be neglected and we can set. 



l-g2-i(l + g2 



(5.6) 



to a very good approximation. 

Using such mode functions we can estimate the quantum fluctuations 5'S(r) for short 
times. Inserting ^Pqir) given by eq. ( |5.6|) in eq.( |3.7| ) the integral is dominated by g = for 
growing times r. 
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Sfr) 



-0 dq 



1 + ^(1 + ?^) 



.2t 



8r3/2 



(5.7) 



Using this estimate for the quantum fluctuations S(r), we can now find the value of the 
spinodal time scale ti at which the back-reaction becomes comparable to the classical terms 
in the differential equations. Such a time is defined by 5fS(ri) ~ 1. From the results 
presented above, we find 



log 
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+ ^ log log 
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V TTfif 




9. 





(5.8) 



The time interval from r = to r ~ ri is when most of the particle production takes 
place. After r ~ ri the quantum fluctuation shut-off the exponential growth of the modes 
and particle production slows down dramatically. 

We list here the values of the spinodal time according to eq.( ^.8|) and the corresponding 
position of the maximum of 5'S(r) for different values of g = 10^" and a fixed r]o = 10~^. 



n = - logio 9 


Spinodal time (ri) 


Numerical Result 




estimate 


for Ti 
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7.85 


7.87 
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9.14 


9.14 


7 


10.40 


10.40 


8 


11.65 


11.40 


9 


12.89 


12.87 


10 


14.14 


14.10 


12 


16.55 


16.90 


14 


19.00 


19.29 



We see that the two values are equal up to 



70. 



B. Asymptotic nonlinear evolution I: numerical analysis 



We present in this section the numerical analysis for the time evolution after the spinodal 
time Ti when the quantum backreaction from 5'S(r) becomes important and a full solution of 



the non-linear equations ( |5.1[ - |5.2| ) is required. We have carried out a numerical integration 
of the integro-differential equations for a wide range of initial amplitudes and couplings 
implementing the same algorithms and with the same precision as in the unbroken symmetry 
case before. In the broken symmetry case we find numerically that the effective mass squared 



(5.9) 



vanishes for r ^ oo. 

A detailed analysis that includes FFT, reveals that Ai^^r) tends to zero oscillating with 
decreasing amplitude. For r > ri a very precise fit is obtained in the form [see fig. 15], 



(r) = - sir) 

T 



COS 



2r 



2c2s(r) log ( — ) + 7 



(5.10) 
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where A, C2 and 7 are constants and s(r) is a slowly decreasing function with s{ti) ^ 1. 
This slowly varying function can be well approximated by: 



.-{r/r2r 



(5.11) 



For small initial amplitudes (r^o < 10 A and tq are independent of rjQ and their dependence 
on g can be numerically fit as follows, 

, ro~l/v^ , A ^ 0.64 log- - 1.9 , C2 ^ 0.6 - 0.16 log- . (5.12) 

9 9 

for small couplings {g < 10~^). 

These expressions reveal that a second dynamical time scale T2 emerges in the broken 
symmetry case. Furthermore, this new scale is widely separated from the non-linear scale, 
i.e. T2 ^ Ti for weak coupling and small initial amplitude. 

We also notice that the frequency of the oscillation of the effective mass is given by 
the initial mass (for very small tjq), which is in agreement with the situation in the un- 
broken phase [see eq. ([4.3|) ] since in this broken symmetry case the effective mass vanishes 
asymptotically. 



Power law relaxation was also reported in [34] for the broken symmetry phase. 



C. The evolution of the mode functions 



Since the effective mass vanishes asymptotically, the g-modes (/'^(t) oscillate as 



(5.13) 



behavior that is confirmed in our numerical analysis. 

As for the unbroken symmetry case, it is convenient to introduce slowly varying ampli- 



tudes Ag{T) and -Bg(T), 



Mr) = ^e 



Bqir) ^ - e 



-tq T 



+iqT 



(^g(r) + - V3g(r) 



(5.14) 



to separate the evolution on short time scales from that on the long time scales. We can 
thus express the mode functions ^qij) in terms of Aqij) and Bqij-) as follows. 



^q{r)=Aq{T) e"'^ + i?,(r) e^"- . 
We find from eqs. (|OB|)-(CT), 

lim A,(r) = A,, 



(5.15) 



(5.16) 



lim 5,(r) = Bq 



(5.17) 
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We obtain from eq.( |5.14|) for the square modulus of the modes, 

|<^,(r)|2 = |A,(r)|2 + |5^(r)|2 + 2|A,(r) 5,(r)| cos [2gr + 0,(r)] (5.18) 
where we set 

A,(r)i?,(rr = |A,(r)i?,(r)|e^^^M. (5.19) 

The wronskian relation ( |3.6|) and eq.( ^.13|) imply that the functions Aq{T) and Bq{T) are 
related through 

\B,{r)\' - \Mr)\' = . (5.20) 

plus terms that vanish fast asymptotically. We plot in fig. 18 the (scaled) modulus as a 
function of time 



M,{t)^^^\Mt)\^ + \B,{t)\^ (5.21) 

and (pqir) for some relevant cases. This figure shows that as anticipated, Mq{T) and (pqir) 
vary slowly with r. 

For small coupling g, |(y9q(r)p, |i?g(r)p and |y4q(r)p are of order 1/g for q in the spinodal 
band < g < 1 and times later than ri p0|j27[l . Therefore, Mg(r) is in that case of order 
one. Moreover, eq.( |5.2U| ) implies that |i?g(r)p = |74q(r)p[l + 0{g)] and we can approximate 
eq.( |5.18| ) as follows, 

g\^q{T)\^ = Mqirf {1 + cos [2gr + 0,(r)]} [1 + 0{g)] , (5.22) 

for < g < 1. 

We see from the numerical results (fig. 16) that the momentum distribution of produced 
particles evolves with time. In particular, the position of its peak g_B(r), decreases with 
time. From the numerical results we can fit the peak position by the form: 

gB(r) ^ , r > n . (5.23) 

T -Tb 

where ri ^ 3.8 for g < 10"'^, r^o < 10~^ and tb ~ Ti. 

Since Al^(oo) = 0, the asymptotic equation for the zero mode, fjij) = 0, only admits a 
constant as bounded solution. This is precisely what the numerical analysis yields. 

We find two non-linear resonant bands for t > ti. The first for < g < ^^(t) and the 



second for qdr) < q < yl — ?7o ~ 1. The value of gc(T) follows by perturbation analysis 
[see next section]: 

qcir) ^ 1 ^ sir) , r > n . (5.24) 

T-Tc 

where rc ~ ri. 

We see from eqs.( |5.23"D and ( |5.24| ) that both unstable bands shrink to zero for r cxd. 
The modes in these nonlinear resonant bands grow as powers of time. The mode exactly 
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at g = is not resonant. ri{T) tends to a constant value for r — oo while ipq=o{T) grows 
linearly with the time. 

The order parameter ri{T) tends for late times to a nonzero limit that depends on the 
initial conditions. 

The modes between the unstable bands [^^(r) < q < q'c(T)] oscillate in time with 
constant amplitude [see figs. 19]. 

Most of the particles are produced just before the time ri (see fig. 17). Then, the particle 
number oscillates with decreasing amplitude around its asymptotic limit. This limiting value 
is proportional to and independent of r/o for small tiq and g and approximately given by: 



Ar(oo) 



0.330 



(5.25) 



In the present case the particle number saturates around ti. 

The vanishing of the effective mass in the r — > oo limit implies that the g-modes become 
effectively free. As for unbroken symmetry [see sec. IIIA], the system tends asymptotically 
to a limit cycle. All modes oscillate harmonically for large enough times, since for t = oo 
both unstable bands (0 < g < ^^('r) and qc{T) < q < 1 shrink to zero). In addition, the 
expectation value tends to a small constant value 77(00) = rjoo- 

In the present broken symmetry case we find that the mode q = grows linearly in time, 
whereas modes with q = \Jl — r]^ 1, grow with a power law for r > ri approximately 
given by 



=1 ^ 



r 



B. 



q=l\ 



-IC2 



ci 



(see fig. 18) with 



ci ^ 0.24 ; C2 ^ 0.6 - 0.16 log 



1 



(5.26) 



(5.27) 



A more precise numerical fit yields that the dependence of ci on g can be bounded by 

Aci 



Alog^- 



< 4. 10^ 



(5.28) 



D. Asymptotic nonlinear evolution II: perturbative and multitime scale analysis 

As in the unbroken symmetry case, we can determine analytically the position of the res- 
onances by performing a perturbative analysis in a similar manner, in terms of the advanced 
Green's function. 

The effective mass 7Vl^(r) vanishes for r — > cxd in the broken symmetry case [see 
eq. (|5.10|) ]. Because the function s(r) varies on even longer time scales of order l/^/g we 
can perform a perturbative expansion on intermediate asymptotic scales, between ri and 
the much longer time scale associated with s(r). In this intermediate asymptotic regime we 
can consider s(r) to be constant. 

We can thus study the mode equations 
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+ q' + M\r) 



Vqir) = 



(5.29) 



for late times considering A4'^{t) as a small perturbation. That is, we can write the integral 
equation including the boundary conditions ( p.l3| ) 



°° , , smgfr 
dr — 



tq T 



(5.30) 



Iterating this integral equation and using eq. (|5.1(]|) yields after calculation if we take the 
function sir) as a constant, 



Ag e*«" + Bg e^*^" 



ir 



Bg e"*"" 



Ag e^" 



q+l 



g-1 



9-1 



(5.31) 



where iIj{t) = 2 r + 2 C2 s(r) log (^;^^ + 0. 

These expressions are singular at g = ±1 revealing the the resonant band below q = 1. 
Actually, these singularities disappear when the decrease of s(r) with r is taken into account 
on time scales r >> T2. 

Notice that there is no resonance at g = for broken symmetry. 

For q'^ ^ 1 the perturbative expansion breaks down for the intermediate asymptotic time 
scales Ti << r << T2. 

In order to understand the behavior of the mode functions in the resonant regions, we 
implement a multitime scale analysis in the same manner as in the unbroken symmetry 
case introducing the parameter e = 1/ri and introducing the long time scale Ti = eTo with 
To = r. The analysis is carried out in exactly the same manner as before. The result is 
that the amplitudes Aq=i ; Bg=i are exactly given by eq. (|5.26| ). Furthermore a consistency 
condition arising from the multitime scale analysis is that 



A = 4^ci + 4^4c2 . 
This condition is well verified in the numerical calculations. 



(5.32) 



Since the effective squared mass vanishes asymptotically and there are no resonances at 
g = for sufficiently late times the q = mode behaves as 

y,o(r)=L + Mr (5.33) 

where L and M are complex coefficients that can only be obtained from the full time 
evolution. The wronskian relation ( p.6|) implies that 

/m[LM*] = l, (5.34) 
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showing that M 7^ 7^ L. 

It follows from eq.( p.8|) that the number of g = quanta grows as for asymptotically 
large r, 



No{r) 



2 _2 



1 + ?7g |Mp r 



(5.35) 



Notice that the total number of particles tends to a constant for t —>■ 00 because Nq{t) 
does not contribute to the total due to the vanishing of the phase space factor as at low 



momentum. This linear growth of the homogeneous quantum mode was also noticed in Mi 



This situation is similar to Bose-Einstein condensation in which case the excess number 
of particles at a fixed temperature goes into the condensate, whereas the total number of 
particles outside the condensate is fixed by the temperature at zero chemical potential. The 
g = mode does not contribute to the energy, the pressure and the total particle number 
but it will become macroscopically occupied for r y/V with V the volume of the system, 
in which case the number of particles in the zero momentum mode becomes of the order 
of the spatial volume. In this case this mode must be isolated and studied separately from 
the q modes because its contribution to the momentum integral will be cancelled by 
the small phase space at small momentum much the same way as in usual Bose-Einstein 
condensation. 

The presence of a Bose condensate through this macroscopic zero momentum mode sig- 
nals spontaneous symmetry breaking even when the order parameter remains zero. Therefore 
we identify the linear growth in time of the g = mode as the onset of a novel form of Bose 
condensation of Goldstone bosons and symmetry breaking in the collisionless regime and 
in the absence of thermalization. This form of Bose condensation of Goldstone bosons in a 
collisionless regime is similar to that reported recently within a different context . 

The time scale for the formation of the condensate 



(r ^ y/V ), is obviously much larger 
than the nonlinear scale ti and the thermalization scale when corrections of order are 
included. Therefore, the formation of the Bose condensate including corrections, which 
will include collisions, will require a further understanding of the dynamics and of the time 
scales involved. 



E. Energy and Pressure for Broken Symmetry: sum rules 



As shown in refs. [^,0 the energy ( 4.63| ) and the pressure ( [4.65 ) can be rewritten for 
broken symmetry as follows 



2\Mn\^ 
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Pir) = g f q^dq 
Jo 



3 



+ r,' + 0{g) - e , (5.36) 

(5.37) 



up to corrections of order ?7g << 1 where qu is the maximum spinodally unstable wavevector 
(for r]l << 1 g„ = 1. 

As for the unbroken symmetry case, we can use the asymptotic value of the effective mass 
squared and the conservation of the energy to compute the first two moments of M^(oo) 
and to establish similar sum rules. 

The vanishing of A^^(oo) implies that gT,{oo) = 1 + C(?7o) and hence 

J\'dq M2(oo) + 0{g,ril) = 1 . (5.38) 



where we used eqs.( |3.7| ) and ( |5.22|) . 



Furthermore, equating the expression for the energy at r = oo [eq. (|5.36| )] with its initial 
value 

e = l + Oiv'o) (5.39) 

yields the second sum rule: 

j\UqM;{^) = ^ + 0{g,r^',). (5.40) 
Using these sum rules allows to compute the pressure in the r — > oo limit from eq.( p.36| ) 



p(oo) = ^ + 0{g,vl) . (5.41) 



leading to the equation of state for radiation 

p(oc) = ^e + 0{g,r]l) (5.42) 
as is expected since the bulk of the produced particles are massless. 

VI. CONCLUSIONS AND FURTHER QUESTIONS: 

In this article we have studied both numerically and analytically the asymptotic non- 
equilibrium dynamics of relaxation in a scalar field theory in the collisionless regime. We 
have focused on the relaxation of initial states of very large energy density that require a 
non-perturbative treatment in a controlled manner, maintaining renormalizability, energy 
conservation and all of the relevant conservation laws. Detailed numerical analysis revealed 
the presence of new dynamical and non-perturbative time scales and non-linear resonant 
bands that result in the power law growth of quantum fluctuations and a power law of 
relaxation for the expectation value of the scalar field. A separation between the time scales 
for weak coupling allowed us to implement a dynamical renormalization group resummation 
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of secular terms via the method of multi-time scales and confirmed and completed the 
numerical results. 

This dynamical renormalization group resummation of secular terms leads to power 
law relaxation with anomalous dynamical exponents which are non-universal depending 
non-perturbatively on the coupling constant. 

In the unbroken symmetry phase the expectation value vanishes asymptotically, transfer- 
ring all of the initial energy into production of massive particles, as a result of the non-linear 
resonances and despite the presence of perturbative thresholds to particle production. The 
asymptotic distribution is non-thermal non-perturbative in the band where parametric res- 
onance takes place at early times. It can be interpreted as a 'semiclassical condensate' in 
the unbroken phase, a result of the relaxation of the initial energy density. The equation 
of state of the produced particles interpolates between dust and radiation domination as a 
function of the initial amplitude of the expectation value of the scalar field. 

In the broken phase the numerical evolution revealed a hierarchy of time scales, and the 
relaxation of the order parameter is with anomalous power laws. Again an implementation 
of the dynamical renormalization group revealed the presence of non-linear resonances that 
result in particle production after the spinodal time. Non-universal power law relaxation 
appears with exponents depending non-perturbatively on the coupling constant. The cou- 
pling constant dependence is similar for the unbroken and broken symmetry cases. The 
asymptotic distribution of particles is localized at low momenta, is non-perturbative and 
non-thermal, but the equation of state is that of radiation. We also found the onset of a 
novel form Bose condensation in the coUisionless regime and without thermalization but 
with extremely long time scales. 

Although this body of results provides for a deeper understanding of the non-equilibrium 
dynamics in the coUisionless regime, and thus we believe it represents a quantum field 
theory example of non-linear dynamics there are very many unanswered questions that 
deserve further study, i) the fact that the asymptotic distributions arc of 'soft' momenta 
and semiclassical in the sense that the amplitude of the particle distribution inside the 
bands is oc 1/A, and that there is a separation of time scales for weak coupling should allow 
a 'coarse grained' description in terms of quasi-particle distributions a la Boltzmann. ii) 
The next order in 1/N must be pursued to incorporate consistently collisional processes, 
probably i) must be understood before this to separate the microscopic from the long time 
scales. Collisions will compete with the coUisionless processes and affect the onset of Bose 
condensation found in the coUisionless regime. Here there is the possibility that whereas 
thermal equilibrium is established on collisional time scales, chemical equilibrium may be 
established on much longer time scales resulting in a non-vanishing chemical potential. 

Only a deeper understanding of i) and ii) will lead to a complete understanding of the 
different regimes: coUisionless, kinetic and hydrodynamic that could account for thermal- 
ization, chemical equilibration and all of the important relaxational phenomena. 
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Figure 3: The total number of produced particles as a function of time for r/o = 1, 
g = 10^^^. After the exponential increase around r = ri = 163.7.., M°''^{t) keeps growing. 
For times r > 200, eq. (f4.30|) gives a very good approximation to the numerical results (after 
averaging over oscillations). 
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Figure 6b: The effective mass squared minus its value at r = oo times r as a function of 
time for rjQ = 2, g = 10^^. Tliis function oscillates in time with constant amplitude Ki + K2 
and frequencies 2Moc = 2w^l + t]^/2 and 2Mo = 2J1 + t^ [see eq.(^3D]. 
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Figure 7a: Momentum distribution of the produced particles at r = 25. for rjo = 4, 
g — 10'^. Notice that the single peak present for times r < ti splits into two due to the 
nonlinear resonances. Many more peaks appear for subsequent times as shown in fig. 7b. 
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Figure 7b: Momentum distribution of the produced particles at r = 200 for tjq = 4, 
g = 10^^. Notice the main peak at g = 0.549 associated with the main non-hnear resonance 
(g = 0) and the secondary peak at g = 2.77 associated to the non-hnear resonance at 
q = r]Q/\/2. The position of both peaks are correctly estimated by eqs.(|4.(j|) and ( [4.43p , 
respectively. 
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Figure 8b: 

The amplitude Mq{T) of the mode function (pg^Q^^r) as a function of time for rjo — 4, 
g = 10~^. This amphtude grows faster for r > ri than the q — amphtude but much slower 
than the exponential growth in the forbidden band < q < rjo/y/2 for r < Ti. 
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Figure 9d: 

The phase ^^(t) of the mode function (fiq=o.6{T) as a function of time for rjo — 4:, g — 10"^. 
This phase becomes an oscillatory function as the same time as the modulus Mq{T) [see fig. 
8d] stops to grow. 
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Figure 11a: 

The modulus Mg(r) of the mode function v'g=2.75(T) as a function of time for rjo = 4, 
g = 10~^. This function grows till q = 2.75 gets out of the second nonlinear resonance band. 
The estimate ( [4. 431 ) (r ~ 140.) is in very good agreement with the numerical results plotted 
here. 
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Figure 15a: 

Modulus of the function Aq{Ti) obtained through multi-time scale analysis as a function 
of q^r/ M.OO ■ We give in eg. ( |4.59|) its expression in terms of Whittaker functions. Notice 
the resemblance with the full numerical solutions displayed in figs. 8c and 8d. 
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Figure 15b: 

Phase of the function Aq{Ti) obtained through multi-time scale analysis as a function of 
g^r/A^oo • We give in eq. (|4.59| ) its expression in terms of Whittaker functions. Notice the 
resemblance with the full numerical solutions displayed in figs. 9c and 9d. 
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Figure 17: 

Momentum distribution of the produced particles at r = 200 for broken symmetryr/o = 
10"^, g = 10^^^. Notice the main peak at g = 0.020 and the secondary peak at g = 0.946 
associated to the non-hnear resonance at g = 1. The position of both peaks are correctly 
estimated by eqs.( |5.23|) and (|5.24| ), respectively. 
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